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Abstract 

We report improved measurements of exclusive two-body charmless hadronic B meson decays 
B — ► n'h, where h is a charged kaon or pion or a K°. These results are obtained from a data sample 
that contains 386 million BB pairs collected at the T(4S) resonance, with the Belle detector at the 
KEKB asymmetric energy e + e~ collider. We measure B(B° -► n'K°) = (56.6+| | ± 3.3) x 1(T 6 , 
B(B + -► n'K + ) = (68.6 ± 2.1 ± 3.6) x 1(T 6 and B(B+ -► rfir + ) = (1.73 ± 0.12) x 10~ 6 , 
where the first and second errors are statistic and systematic, respectively. The CP asymmetries 
in the charged modes are measured and no evidence for direct CP violation is found. We measure 
A C p{B ± -► jfK*) = 0.03 ± 0.03 ± 0.02 and A CP {B ± -► t/tt*) = 0.15±g;|| 

PACS numbers: 
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The properties of the decay B — > are still to be understood by theory. The channel 
B — > v(K has the largest branching fraction of all charmless hadronic B decay modes. In 
the Standard Model (SM]the decay B — »■ rj'K is thought to proceed dominantly via gluonic 
penguin processes The diagrams |4| that contribute are shown in Fig. ^ The 

measured branching fractions for B — > rj'K 5 J(| 0| are larger than expectations from the 
generalized factorization approach 0, 0, [hJ 1 1| . This has led to speculations that SU(3)- 
singlet couplings unique to the rj' meson or new physics Q, E3, H contribute to the 
amplitude. 




u. d 



K+,K° 



a . d 



U . d 



a . d 



B+,B® 




u c t 



u, d 



K + ,K a 



(a) 



(b) 




FIG. 1: Feynman diagrams describing the decays B — > rj'K + '°: (a), (b) internal penguins, (c) 
external tree, (d) flavour-singlet penguin. 

The B + — > 7]'tt + decay proceeds via similar processes as B — > rj'K with b —> u(d) tree or 
penguin decays. However, for the latter all contributions are additive, while for the former 
tree and penguin contributions interfere destructively. 

In the SM, direct CP violation arises from the interference of two or more amplitudes 
with different strong and weak phases j^j. Many charmless hadronic B meson decays 
contain both tree and penguin amplitudes and provide a rich sample for direct CP violation 
studies. The decay B — > tj'tt is expected to have comparable contributions from both tree 
and penguin diagrams, and therefore may have a large asymmetry. For the decay processes 
dominated by one single amplitude in the SM, a non-zero CP asymmetry may indicate 
additional amplitudes in the decay and hence provide a hint of new physics. This is also the 
case for B — ► rj'K, which is thought to proceed entirely through the penguin amplitude and 
no CP asymmetry is expected. 
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With the convention that a B meson contains a b quark, the direct CP asymmetry in a 
B — > / decay is defined as 

r(g -> /) - r(g -> /) 

CP T(B^f)+T(B-+ /)' 

where / is the final state and / its CP conjugate. In the experiment we measure the 
branching fraction BiB — > f) and B(B — > /), which are proportional to the partial widths. 

In this paper, we update the measurements of B — > rj'K + ,rj'K° ,r]'7r + with a sample 35 
times larger than our previous dataset We also report the charge asymmetry of self- 
tagged decays. All results are based on a data sample that contains 386 million BB pairs 
collected with the Belle detector at the KEKB asymmetric-energy e + e" (3.5 on 8 GeV) 
collider jH]. KEKB operates at the T(4S f ) resonance (yfs = 10.58 GeV) with a peak 
luminosity that exceeds 1.5 x 10 34 cm~ 2 s -1 . 

The Belle detector is a large-solid-angle magnetic spectrometer that consists of a silicon 
vertex detector (SVD), a 50-layer central drift chamber (CDC), an array of aerogel threshold 
Cerenkov counters (ACC), a barrel-like arrangement of time-of-flight scintillation counters 
(TOF), and an electromagnetic calorimeter (ECL) comprised of CsI(Tl) crystals located 
inside a super-conducting solenoid coil that provides a 1.5 T magnetic field. An iron flux- 
return located outside of the coil is instrumented to detect K\ mesons and to identify muons 
(KLM). The detector is described in detail elsewhere j22|. Two inner detector configurations 
were used. A 2.0 cm beampipe and a 3-layer silicon vertex detector was used for the first 
sample of 152 million BB pairs (Set I), while a 1.5 cm beampipe, a 4-layer silicon detector 
and a small-cell inner drift chamber were used to record the remaining 234 million BB pairs 
(Set II) 0. 

The rf mesons are reconstructed via either rf — > 777r + 7r~, with rj —>■ 77, or rf —>■ p°7. K° 
mesons are reconstructed only via K$ — > 7r + 7r~. 

Primary charged tracks were selected with dr < 0.5 cm and \dz\ < 2.0 cm, where dr and 
dz are the impact parameters perpendicular to and along the beam axis, with respect to the 
run dependent interaction point (IP). 

Likelihoods for kaon and pion hypotheses, Ck and C n , are obtained by combining in- 
formation from the CDC (dE/dx), ACC and TOF systems. The likelihood ratio TZk = 
£k/(£tt + £-k) ranges between (pion-like) and 1 (kaon-like) and we require TZk > 0.6 for 
kaons and TZk < 0.4 for pions from Kg, which keeps about 86% of the kaon/pion candi- 
dates, and TZk < 0.9 for pion candidates from the 7/ decay with an efficiency of about 98%. 
Prompt pions in B + — ► r/7r + are required to pass the tight selection, TZk < 0.1, which has 
an efficiency of about 85% and a kaon fake rate of about 4%. 

The 77 meson is reconstructed from two photons, each with an energy of at least 50 MeV. 
The rj mass window is chosen to be 0.5 GeV/c 2 < M(77) < 0.57 GeV/c 2 , which corresponds 
to +2/ — 3 standard deviations (a) from the nominal 77 mass given by Particle Data Group 
(PDG) 0. 

We retain p° candidates in the mass range 550 MeV/c 2 < M(tt + tt~) < 870 MeV/c 2 , where 
M(tt + tt~) is the p° candidate mass. We require the transverse momenta of the daughter pions 
to satisfy the requirement, Pt( 71 ") > 0.2 GeV/c. This suppresses around 40% of background 
while retaining 86% of the signal. 

The requirements for the rj' candidate mass depend on the decay channel. For B — ► rj'K 
decays, we select rf candidates within ±3.4<r window for rj' — > 777r + 7r~ (0.945 GeV/c 2 < 
M(777T7r) < 0.97 GeV/c 2 ) and ±3<r for rj' -»■ p°7 (0.935 GeV/c 2 < M(py) < 0.975 
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GeV/c 2 ). For B + — > r]'7i + decay, the 77' mass windows are tightened to ±2. 5a (0.95 
GeV/c 2 < M^nn) < 0.965 GeV/c 2 and 0.941 GeV/c 2 < M(pj) < 0.97 GeV/c 2 ). A 
weak requirement on the 77 meson decay angular distribution, h(i]) < 0.97, has been applied 
for the rf — > 7]7T + 7r~ mode to suppress combinatorial background, where h{rj) is the cosine 
of the angle between the rf momentum and the direction of one of the decay photons in the 
r) rest frame. 

Kg candidates are reconstructed from a pair of oppositely charged particles with invariant 
mass within the range 485 MeV/c 2 < M(7r + 7r~) < 510 MeV/c 2 . We require the vertex of a 
K s to be well reconstructed and displaced from the interaction point and the K s momentum 
direction to be consistent with the K s flight direction. 

B meson candidates are then reconstructed combining an rf meson and one of the h 
candidates. Two kinematic variables are used to extract the B meson signal: the energy 
difference AE = Eb — -Ebeam and the beam-energy constrained mass Mb c = \J E^ e&ul — P§, 
where i?beam is the beam energy and Eb and Pb are the reconstructed energy and momentum 
of the B candidate in the T(AS) rest frame. The events that satisfy the requirements, 
Mb c > 5.2 GeV/c 2 and \AE\ < 0.25 GeV are selected for further analysis. 

For events with multiple B candidates, the best candidate is selected based on the \ 2 
of a vertex fit, that optimizes a vertex for all charged tracks in the final state, and a mass 
X 2 - The latter term is necessary because of the photons in the rj', which may introduce 
additional multiple candidates even for the same set of charged particles. The formula used 
to calculate the x 2 i s: 

X 2 = Xvertcx + 1(M(V') - 77v)/<ty] 2 , (1) 

with Xvertex being the \ 2 from the charged particle vertex fit, M(rj') the rf candidate mass 
and mrf the nominal mass of the rf and cy = 0.008 GeV/c 2 the width of the rf mass 
distribution. About 10% of events have multiple candidates and of these 10% are due to 
multiple photons. 

Several event shape variables (defined in the center of mass frame) are used to distinguish 
the more spherical BB topology from the jet-like qq continuum events. The thrust angle 
9t is defined as the angle between rf momentum direction and the thrust axis formed by 
all tracks not from the same B meson. Jet-like events tend to peak near | cos9t\ — 1, while 
spherical events have a flat distribution. The requirement | cos#r| < 0.9 is applied prior to 
any other event topology selections. 

Additional continuum suppression is obtained by using modified Fox-Wolfram mo- 
ments (HI and the angle 9b between the flight direction of the reconstructed B° candidate 
and the beam axis. A Fisher discriminant (J 7 ) [26] is formed by a linear combination of 
cos#t, S± and five modified Fox- Wolfram moments. S± is the ratio of the scalar sum of the 
transverse momenta of all tracks outside a 45° cone around the rf direction to the scalar sum 
of their total momenta. Probability density functions (PDFs) are obtained from signal and 
background MC data samples. These variables are then combined to form a event topology 
likelihood function C c =P c (cos9b) ■ Pci^ 7 ), where P c = PDF of signal (s) or continuum 
background (qq). Signal follows a 1 — cos 2 9b distribution while continuum background is 
uniformly distributed in cos9b- We select signal-like events by requiring a likelihood ratio 
IZc = C s / (C s + C q q) criteria optimized by MC studies to suppress continuum background. 
For channels with an rj' — > p°7 decay an additional variable cos 9-h, which is the angle be- 
tween the r/' momentum and the direction of one of the decay pions in the p rest frame, is 
included for better signal-background separation. 
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Further background discrimination is provided by the quality of the B flavor tagging of the 
accompanying B meson. We use the standard Belle B tagging package j2?| , which gives the 
B flavor and a tagging quality r ranging from zero for no flavor and unity for unambiguous 
flavor assignment. We divide the data into three r regions and separately optimize the TZc 
requirements using signal and continuum background Monte Carlo samples. 

The signal yields (N$) are extracted using extended unbinned maximum-likelihood fits 
to two-dimensional (AE,M^, C ) distributions. An extended likelihood function is: 



L(N s ,N Bj ) = 

— TT 

N\ M 



N s P s (AEi,M hCi ) + ^2N B .P Bj (AE i ,M, 



be, 



(2) 



where iV is the total number of events, i is an index running over the events and Ps and 
P Bj are the PDFs for signal and background, respectively, and the index j runs over all 
background sources. The signal yield N s and background contents N B . are determined 
by maximizing the L(N$, N Bj ) function, where the variable N B . defines a j-dimensional 
submanifold of all different backgrounds. 

To take into account the efficiency differences in two data sets, Ns is calculated as: 

Ns = 61 N m B + e 2 N BB - 2 B, (3) 

where B is a branching fraction, and the and N BB -_ are the efficiency and the number of 
BB pairs for Set I and Set II. We assume that the numbers of B + B~ and B°B° pairs are 
equal. 

For charged B decays, we divide the data into two samples for positive and negative 
charges and extract the charge asymmetry Aqp in addition to the branching fraction B, 
using the formula: 

N ± = 0.5 (It A C p)N s , (4) 

with N + (iV_) the number of positively (negatively) charged kaons or pions and Ns the 
number of signal events. 

The PDF shapes for each contribution are determined by MC studies. We assume the 
signal shapes for AE and Mb c to be inde pen dent and model the signal using a Gaussian 
with an exponential tail (Crystal Ball line) [28[ plus a Gaussian for AE and a Gaussian with 
an exponential tail for Mb c . The shape parameters are fixed from the signal MC sample, 
except for the means and widths of both AE and Mb c distributions, which are determined 
in the fit of B — > rj'K for decays with K + or K s in the final state for each data set (Set I 
and Set II) independently. For B + — > i]'7i + we use the parameters obtained from the fit to 
the charged kaon mode. 

We consider up to four types of backgrounds separately in the fit: continuum, b —> c, two 
types of charmless decays mentioned below. The AE and M^ c distributions for continuum 
backgrounds are found to be largely uncorrelated and are thus modeled with two independent 
one-dimensional functions. For AE we assign a first or second order polynomial and the 
Argus function [29] is used to model the M bc distribution. Charmless B decays and b — ► c 
backgrounds are modeled with 2-dimensional smoothed histograms. The contributions from 
charmless B decays are modeled with two smoothed histograms, one for one decay that 
makes a large contribution to the background, and one for all other charmless decays. For 
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B + — > f]'K + , the dominant decay decay mode that is modeled separately is B — > t]'K*; for 
5° -> 7/ifg it is B -> p°Kg; and for 5+ -> 7/71-+ it is the 73+ -> t/^ + feeddown. The 
feeddown in B + — > 7/71- + is modeled with the same PDFs as used for the signal, shifted in 
AE and with an additional correction factor for the change in the width in AE. 

The number of signal and background events are free parameters in the maximum likeli- 
hood fit. The signal mean and width shape parameters for kaonic decays and the continuum 
shape parameters are also free in all fits. For 7/ — > 777r + 7r~ modes our background MC 
studies show that no contributions from b —>■ c decays are expected. For B —>■ rj'K charm- 
less backgrounds are also not expected. For B + — > i]'tt + there is a small contribution from 
charmless B decays in addition to the B + — ► i]'K + feeddown. For 7/ — > p°7 modes, we float 
the number of events for b —>■ c and two charmless B contributions. Other parameters are 
fixed to values determined from MC or sideband data studies. For the combined results, we 
use a simultaneous fit with the branching fraction and the charge asymmetry as common 
parameters for rf — > 777r + 7r~ and 7/ — > p 0/ y. The projection plots of the fits are shown in 
Fig. 120 The reconstruction efficiencies for each decay and results of the fits are displayed 
in Table [I] The upper part of the table displays fit results for individual fits for the two 
subdecay modes 77' — > rjir + ir~ and 7/ — > p°7. In addition, we list the charge asymmetry of the 
continuum background obtained from the fits. Since continuum background is thought to 
be distributed equally for both charges, any significant observed asymmetry would indicate 
a bias and is therefore included in the systematic error (see below). 

We find the significance of B + — > 7]'tt + yield is 3.0, which is calculated as a = 
^2 ln(L max /L ), where L max and L denote the maximum likelihood value and the likelihood 
value at zero branching fraction, respectively. The systematic error (mentioned below) is in- 
cluded in the significance calculation by substracting the systematic error from the obtained 
branching fraction and recalculating the significance. 

The reconstruction efficiencies are determined from signal MC samples before including 
subdecay branching fractions and are around 16-25% for decays with rj' — > 777r + 7r~ and 9-12% 
for 77' — > p°7 before considering subdecay branching fractions. The efficiencies are calculated 
separately for both Set I and Set II. The efficiency for Set II is typically about 0.5% larger 
than for Set II. Correction factors due to differences between data and MC are included for 
the charged track identification, photon, tt° and rj reconstruction, resulting in a correction 
factor of ~ 0.9. The corrections were determined from detailed studies that are discussed in 
the section on the systematic error below. 

We calculate a goodness of fit (gof) based on binned projections of the 2-dim. fit into 
AE. We define: 

gof = X 2 /dof, (5) 

with x 2 = E^^^ the sum over all bins of the projections with and U{ the number 
of expected and observed events in each bin, and dof the degrees of freedom. Since we use 
an unbinned simultaneous fit in AE and Mb c the projected gof value is not a rigourous 
estimator but can still be used to qualitatively evaluate how well we fit the data. 

Systematic errors are estimated with various high statistics data samples. The sources 
and their contributions are listed in Table ITJl The dominant sources of background are 
the uncertainties of the reconstruction efficiency of charged tracks, the uncertainties in the 
reconstruction efficiencies for 77 mesons and photons and the uncertainty of the PDFs shape 
and other parameters, which are estimated by varying each parameter of the PDFs by the 
la uncertainty in its nominal value. The changes in the yield are added in quadrature. 
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TABLE I: Signal efficiencies with (etot) and without (e) subdecay branching fractions included and 
averaged for Set I and Set II for rj — > rjir + ir~ and rj — > p°7, branching ratios B, asymmetry Aqp 
for signal and background and goodness of fit (gof). The errors in e are statistical uncertainties 
in the MC while those in etot include errors of secondary branching fractions. For others, the first 
errors are statistical and the second (if given) are systematic errors. 





B^ — > rj 


B v — ► 77 il 




e[r] — > 7/7T 7T ) [7o\ 


ZO.ZO ± U.ll 


iy.y ± U.ib 


Lb. 6 ± U.iU 


, /„/ . „_+_-n ran 

Ctat{V — > rjTT 7T ) [/o\ 


4.41 ± U.14 


l.ZU ± U.U4 


9 4- n 1 n 

Z.oO ± U.IU 


yieiu {ij —> i/Ji n ) 




04Q o i oi q 

Zj l ±*J.Zj _1_ ^r±.t7 


-L 1 . t: _1_ 1 .W 


D{T] —> T]7T 7T ) [1U J 


u / .1 zc z.u 


52 5+ 4 - 8 


l.OO_ 70 




—0 003 + 036 






e („' _» n0 7 ) [%1 


10.06 ± 0.09 


11.6 ±0.19 


9.5 ±0.08 


6tot(V- A) [%] 


2.97 ±0.10 


1.18 ±0.04 


2.80 ±01.0 


yield (rj' — » p°7) 


823.3 ±42.9 


278.3 ±24.1 


23.9 ±10.1 


W -> A) [10- 6 ] 


71.9±|;f 


6i.HI:! 


2.21 ± 1.28 


A C p(r?' -► A) 


0.09 ±0.05 




—0 D9 +0 - 66 
u.uy_ 74 


yield 


1952.2 ±60.6 


519.9 ±32.3 


37.8 ± 14.5 


i3[10" 6 ] 


68.6 ±2.1 ±3.6 


0D - D -3. 5-3.2 


l-73±g;|| ±0.12 


A C p 


0.029 ± 0.028 ± 0.02 




n 1 c+0.39+0.02 
u - io -0.38-0.06 


Aq p (continuum) 


-0 01 q+ - 008 




o.oi8±g:8§s 


gof 


1.14 


1.51 


1.16 




AE [GeV] M bc [GeV] 



FIG. 2: a) AE 1 and b) M^ c distributions for B + — > rj'K + (for 77' — ► ?7vr + 7r — and 77' — > p°7 
combined) for the region M bc > 5.27 GeV and -0.1 GeV < < 0.06 GeV, respectively. The 
histograms represent data, the red small dashed line the signal contribution, the light blue dotted 
line continuum background, the dark blue large dashed line b — ► c backgrounds and the pink dash- 
dotted line charmless B contributions. The green solid line is the sum of all contributions to the 
fit. 
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AE [GeV] M bc [GeV] 



FIG. 3: a) AE and b) M bc distributions for B° -> r/ifg (for 7/ — > r]ir + ir and 77' — > /9°7 combined) 
for the region M hc > 5.27 GeV and —0.1 GeV < AE < 0.06 GeV, respectively. The histograms 
represent data, the red small dashed line the signal contribution, the light blue dotted line con- 
tinuum background, the dark blue large dashed line b — > c backgrounds and the pink dash-dotted 
line charmless B contributions. The green solid line is the sum of all contributions to the fit. 




-0.25 -0.2 -0.15 -0.1 -0.05 0.05 0.1 0.15 0.2 0.25 



AE [GeV] 



120 



100 




5.22 5.24 5.26 5.28 

M bc [GeV] 



5.3 



FIG. 4: a) AE and b) M bc distributions for B + — > t]'it + (for if — ► r/7r + 7r _ and rj 
combined) for the region M bc > 5.27 GeV and -0.1 GeV < AE < 0.06 GeV, respectively. The 
histograms represent data, the red small dashed line the signal contribution, the light blue dotted 
line continuum background, the dark blue large dashed line b — > c backgrounds, the pink dash- 
dotted line charmless B contributions and the black medium dashed line B + — ► r]'K + feeddown. 
The green solid line is the sum of all contributions to the fit. 



The normalization of the K feeddown to B + — ► i]'tt + is estimated by varying the assumed 
branching fraction for the feeddown by the ±l<r uncertainty in the present measurement 
of B + — > i]'K + , the error arising from differences in the fake efficiency is included in the 
systematic error for the PDF parameters. Corrections for AE and M bc differences between 
data and MC for the feeddown shape of rj'K to B + — ► 7//T + are considered by varying 
the corrections by one standard deviation and refitting. Systematics arising from the TZc 
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selection are studied by varying the TZc selection and by a large B — > Dtt sample, we 
conservatively use the larger error. The uncertainty of the subdecay branching fractions 
is given in the PDG. The number of BB mesons produced at Belle is estimated from the 
number of continuum-subtracted hadronic events. The uncertainty in this value is included 
in the systematic error. The uncertainty of the particle identification is estimated with 
D* + — > D°7i + decays. Other efficiency systematic errors are found to add up to less than 1% 
and are included with a 1% error. All contributions are added in quadrature and we find the 
systematic errors for the three decays to lie between 5% and 7%. For the charge asymmetry, 
efficiency based systematic errors cancel out and therefore we use the high statistics mode 
B + — > 7]'K + to calculate systematic errors and use the same absolute numbers for B — > rj'n. 
The dominant contribution to the Acp systematic error is from the possible difference in 
the detector response for positive and negative charged particles. We estimate this from 
the continuum asymmetry to be 0.02. Other contributions from fitting and normalization 
together result in a systematic error of 0.002. Thus the total systematic error for the Aqp 
measurement is 0.02. For B + — > i]'tt + we assign an additional source of error by adding an 
asymmetry to the B + — > i]'K + feeddown contribution. We find an error of —0.06, which is 
added in quadrature to the other errors. 

TABLE II: A breakdown of systematic uncertainties in percentage for branching fraction measure- 
ments. Systematic errors from tracking, e 7j7r o ;7) , the combined systematics from these three particles, 
particle ID, e, the systematics from efficiency calculations, TZc selection, PDF, from varying the fit 
parameters by one a, AE/ M^ c corrections, which are varied by ±1<j, the B + — > rj'K + feeddown to 
B + — ► r/ir + , B s , the subdecay branching fraction errors and N B jj, the uncertainty of the number 
of BB, are listed. 



Source 


AB/B [%] 
B+ -f r/K+ B + -► rfir + 


B° V 'K° S 


tracking 


3 


3 


4 




3 


3 


3 


part. ID 


0.7 


0.9 


0.6 


e 


1 


1 


1 


n c 


2 


+1.9 
-2.7 


2 


PDF 


+0.5 
-0.7 


+2.6 
-2.0 


+1.2 
-1.3 


AE/ M bc 




+2.1 
-1.9 




r]'K + feeddown 




+3.5 
-2.9 




B s 


1.5 


1.5 


1.5 


N BB 


1 


1 


1 


Total 


5.2 


+7.1 
-6.8 


5.9 



In summary, we report improved measurements with 35 times more statistics of the 
charged and neutral B — > rj'K decay. We find 3a evidence for B + — > r]'iT + and report the 
charge asymmetry for the decay modes B + — > i]'K + and B + — > i]'tt + . The central values 
of our branching fraction measurements are below current PDG values, but are consistent 
within statistical errors. 
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